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Abstract 



The doorway states under consideration are eigenstates of the hamiltonian which 
is the sum of the kinetic energy and the infinite energy limit of the single-particle 
mass operator. Only Hartree diagrams with the free-space nucleon-nucleon forces 
contribute to this limit, and therefore the observed doorway state energies carry an 
important information about both the nuclear structure and the free-space nucleon- 
nucleon interaction. 

1 Introduction 

The experimental data on the quasielastic knockout reactions (p, 2p), (p,pn), (e, e'p) 
etc. leading to the strongly bound hole states of complex nuclei carry an important 
information about both the nuclear structure and the free-space nucleon-nucleon 
forces. Of course such information is contained in all nuclear data, but these ones 
are distinguished by the fact that the above information is obtained by simple means 
thus being highly reliable. The reasons are as follows. 

1. As shown by M.Baranger [|I] the doorway states for the one-nucleon transfer 
reactions are eigenstates of nucleon in the static nuclear field, see Sect. 2, thus being 
solutions of the problem of particle motion in the central potential well. This is one 
of the most simple problems of quantum mechanics. 

2. As a consequence of contemporary ideas about the AW interaction mecha- 
nism, see Sect. 3, the only contribution to the static field of nucleus is provided by 
the Hartree diagrams with the free-space nucleon-nucleon forces: the two-particle, 
Fig. la, three-particle, Fig. lb, four-particle, Fig.lc, etc. 

* E-mail: birbrair@thd.pnpi.spb.ru 
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Figure 1: Hartree diagrams for the static field of nucleus. 

3. The two-particle contribution of Fig. la is the convolution of the free-space 
two-particle NN interaction with the nucleon density distribution in nucleus, and 
therefore it can be determined from experiment. Indeed, the two-particle forces are 
determined by the properties of deuteron and elastic NN scattering phase shifts 
below the pion production threshold, whereas the nucleon density distributions are 
deduced from the combined analysis of the electron-nucleus || and proton-nucleus 
elastic scattering data ||. 

The information about the many-particle contributions to the static nuclear field 
(hence, about the free-space many-particle NN forces) can be obtained by compar- 
ing the observed doorway state energies with the calculations including the two- 
particle contribution only. In this way we found that the free-space many-particle 
interaction includes at least the three-particle repulsion and four-particle attraction, 
see Sect. 4. 



2 The Baranger theorem 

Evolution of the state arising from a sudden creation of particle or hole in the ground 
state of nucleus A is described by the single-particle propagator || 

S(x,x';r) = -t(A \Tij(x,r)^ + (x\0)\A ) = 
= ie(-r)J2^j^j(^ iEjT -tO(r)J2^k(x)^t(x')e' iEkT (1) 

j k 

with 

m 3 {x) = ((A - l)Mx)\Ao) , Ej = £ (A) - £ 3 {A - 1) 

^ k (x) = (A \ij(x)\(A + l) k ) , E k = £ k (A+l)-£o(A) . (2) 

So the propagator describes the evolution of the hole (particle) state at negative 
(positive) t values. According to Eq.(2) the excitation energy region for the A — 1 
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nucleus is 

- oo < Ej < S (A) - £ g (A - 1) (3) 

(S g (A — 1) and £q(A) are the ground-state energies of A — 1 and A nuclei), whereas 
that for the A + 1 nucleus is 

£ g (A + 1) - £o(A) < E k < oo . (4) 

Such energy scale is convenient for us because the two regions do not overlap in 
stable nuclei. 

The Fourier transform of the propagator 



-t-oo 

G{x,x';e) = J S(x,x';T)e ier dT 



y e-Ej-iS V £-E k + i$ ' 1 ' 

(which is referred to as the single-particle Green function) obeys the Dyson equation 

eG(x,x';e) = 5(x — x') + k x G(x, x'; e) + / M(x, xi;e)G(xi,x';e)dxi , (6) 



k x is the kinetic energy and M(x,x';e) is the mass operator. The latter has the 
following general form 

M(x,x';s) = U(x, x') + J2( x ' x 'i £ ) > ( 7 ) 

where the energy-independent part U (x, x') is the static field of nucleus, and the 
energy- dependent one XX X 5 x ''-> £ ) ls responsible for all kinds of the correlation effects 
(Pauli, particle-particle, particle-hole, ground-state, long-rangle, short-range etc.). 
It has the following high-energy asymptotics || [| 

Y,(x,x';e) = + •■■ (8) 

(the dots in the rhs denote the higher-power terms in respect of e^ 1 )- As a result 
the static nuclear field is the infinite energy limit of the mass operator, 

Uix.x') = lim M(x,x':e) (9) 

e— >oo 

the decomposition (7) thus being unambiguous. 

Now let us introduce the single-particle Hamiltonian 

H sp (x,x') = kJ(x-x') + U(x,x') (10) 
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and its eigenstates 

= J n sp (x,x')ipxW)dx' , (11) 

which are those of nucleon in the static field of nucleus. They are not directly 
observed because they are described by only a part of the total nuclear Hamiltonian. 
Their physical meaning is however understood on the basis of the Heisenberg relation 
according to which the infinite e value is equivalent to the infinitely short time 
interval r. Hence, the eigenstates of H sp , Eq.(lO), describe the very beginning of 
the evolution process under consideration thus being the doorway states for the 
one-nucleon transfer reactions. 

To demonstrate this more explicitly let us use the high-energy asymptotics of 
the Green function (5): 

G( x , x >;e) c ^ = b^ + h^l + h^fL + ... , (12) 

£ £ £ 6 

where 

I (x,x') =J2^A^t^) + J2^k(xM(x') (13) 

j k 

h(x,x') = 52E j * j (x)*}{af) + 52E k * k (x)*t{J) (14) 

j k 

i 2 (x,x') =E^^)^ + (^) + E^ 2 ^(^)^(^). (is) 

j k 

As follows from the spectral representation of the propagator, Eq.(l), 

I (x,x') = i(s(x, x'; +0) - S(x, x'; -0)) (13a) 

I^x.x) = - (s(x, x'; +0) - S(x, x'; -0)) (14a) 

I 2 (x,x') = -i(s(x, x'; +0) - S(x, x; -0)) (15a) 

the above sums thus describing the beginning of the evolution process (S— 

— f^f)- Using the definition (10) and the asymptotics (8) the Dyson equation 
(6) may be written in the form 

eG(x, x- £) = 5{x -x') + J (n sp (x, an) + H^il + G ^ x >. £ ) dxi . ( i6 ) 

Putting Eq.(12) into Eq.(16) and equating the coefficients at the same powers of £~ l 
we get 

Y^^Ax^ix') +J2^k(x)^tW) = S(x-x') (136) 

j k 
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^EjVjix^fix') + J2E k *k(x)Vt(x') = H sv {x,x') (146) 
j k 

'£E%9 j (x)*+W +J2El* k (x)*+(x') = H%(x,x')+U(x,x') . (15b) 

j k 

As follows from Eqs. (14a) and (14b) 

H sp {x,x') = - [S{x, x'; +0) - S(x, x'; -0)) . (17) 

So the evolution of the hole (particle) state begins with the formation of the nucleon 
eigenstates in the static field of nucleus, the Baranger theorem thus being proved. 

Now let us discuss the determination of the doorway state energies E\, Eq.(ll), 
from the experimental data. The weights of the doorway component in the actual 
nuclear states are 



2 

(A) - / Va (x)V j>k (x)dx 



S j,k 



Multiplying Eqs. (13b)-(15b) by ip~^(x)ip\(x') and integrating over x and x' (x 
denotes the totality of space and spin variables) we get 

E*f + E4 A) = i (i3c) 

j k 

E^4 A) +E^4 A) = ^ (i4c) 

j k 

E^4 A) +E^4 A) = el + al (15c) 

3 k 

a \ = J i't (x)H(x, x')ip\(x')dxdx' . (19) 

It is remarkable that in contrast to the widths of the Landau-Migdal quasiparticles 
H the dispersion <j\, Eq.(19), depends upon the wave function ip\( x ) rather than 
the energy E\, thus being roughly the same for all doorway states. In such situation 
it is reasonable to identify a with the largest observed width value. The latter is the 
widths of the peaks in the cross sections of quasi-elastic knockout reactions (p, 2p) 
and (p,pn) ]7| |8|] leading to the lsi/2 hole states. According to the above references 
it is about 20 MeV in all nuclei. 

As seen from Eq.(14c) the doorway state energies E\ are expressed through the 
energies and s-factors of the actual nuclear states. In general case the latter ones 
belong to both the A — 1 and A + 1 nuclei, and therefore the s-factors from two 
different reactions, pick up and stripping, are required. The absolute values of the 
s-factors are, however, measured with a rather low accuracy because of both the 
experimental and theoretical ambiguities. For this reason the energies of weakly 
bound states with \e\\ < a are yet unknown (one should bear in mind that the 
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low-lying states of A 1 nuclei are Landau-Migdal quasiparticles || rather than 
the states of nucleon in static nuclear field). 

The situation is more favourable for the states with \s\\ > a. In this case, see Eqs. 
(3) and (4), the actual nuclear states, over which the doorway ones are distributed, 
belong mainly to either the A — 1 nucleus or the A + 1 one, only one term in the 
lhs (the first for hole states and the second for particle ones) of Eqs.(13c)-(15c) 
thus being active. This is just the case for the strongly bound hole states which 
are excited in the quasielastic knockout reactions (p, 2p) and (p,pn) f7|, U For this 
reason the average energies of the peaks in the cross sections may be identified with 
the doorway state energies within the experimental accuracy of 2-3 MeV. We use 
the facts that the cross section of the quasielastic knockout reaction leading to the 
fixed nuclear state is proportional to the s-factor of this state, and the absolute 
values of the s-factors are unnecessary when all states, over which the doorway 
one is distributed, belong to the same nucleus (in this case the relative values are 
sufficient). 

The experimental data of Refs.f?], |J, which are used in the present work, are not 
free of the following possible ambiguity: the energy of the knocked-out nucleon is 
only about 100 MeV in the experiments. This may be insufficient to neglect the final- 
state inelastic interactions leading to an additional excitation of the final nucleus. 
As a result of such excitations the average energies of the peaks may be shifted 
from the doorway ones because the reaction mechanism is not a pure quasielastic 
knockout in this case. For a greater confidence the additional quasielastic knockout 
experiments (p, p'N) or (e, e'N) are desired, in which the energy of the knock-out 
nucleon would be of order of 0.5-1 GeV. We hope that our work will stimulate such 
experiments. 



3 The static field of nucleus 

Consider the high-energy asymptotics of the Feynman diagrams constituting the 
mass operator. Let us begin with those of first order with respect to the free-space 
NN interaction. The Hartree diagrams of Fig.l are obviously energy- independent. 
But this is not the case for the corresponding Fock diagrams resulting from the 
two-particle, Fig. 2a, three-particle, Fig. 2b, four-particle forces, Fig. 2c, etc. Indeed, 
according to the contemporary ideas the NN interaction proceeds via the exchange 
by either mesons in the Yukawa-like models (OBE M, Paris [TIJ, Bonn [|TTJ, OSBEP 



|i~2"lj ) or quarks and gluons in more sophisticated ones. In any case the interaction 



includes both the momentum and the energy transfer. As a result of the latter the 
Fock diagrams have the e~ l asymptotics. Let us demonstrate this for the diagram 
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(a) 



(b) 



(c) 



Figure 2: First-order exchange contributions to the mass operator. 



of Fig.2a, 



M F (x,x';s) = J 



(20) 



using the Bonn B potential [1 1] for the two-particle NN forces. It is the sum of the 
terms 



2a 



Af + q 2 - to 2 J jj 2 + q 2 -uo 2 



i = Tr,r), p,uj,al,a0,5 (21) 



in the four-momentum space, the form of the meson-nucleon vertices and the sign 
being specified by the Lorentz symmetry of the mesons. Both the sign and the 
Lorentz structure are disregarded here because they are irrelevant for the energy 
dependence. Confining ourselves by the monopole formfactor, a = 1, we get 



M F (x,x';e) = g 2 I ^-W-O 



1 



x 



X 



(27T)' 



2^ M (g) 



X 



* k (x)*tW) 



e-Ej + u^q) Y e - E k - u^q) 



1 - (A 2 - ^ 2 ) 



d 

dX 2 



Y £ - Ej + cu A (q) ^ e - E k - cu A (q) 



2u A (q) 

= v^. Mq) = \/ A2 + 1 2 ■ 

In the e — > oo limit this gives 
g 2 8(x — a/) 



(22) 



Mp(x, x'\ e) 



4vr 2 



A 2 - n 2 

uJil) ~ w A (g) ^ A 3 (g) 



. (23) 
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(a) (b) (c) 

Figure 3: Some second-order diagrams for the mass operator. 

The second-order diagrams of Fig. 3 as well as the higher-order ones contain 
the propagators of intermediate states, and therefore they all have at least the e^ 1 
asymptotics. So the only contribution to the nuclear static field, Eq.(9), is provided 
by the Hartree diagrams. 

The two-particle contribution to the static field of nucleus, Fig. la, is calculated 
with two different models for the free-space two-particle NN interaction, both being 
of clear physical meaning and containing small number of adjustable parameters. 
The first, the Bonn [IT|, is the sum of the OBE potentials with the vertex form- 



factors, Eq.(21). The parameters are adjusted to reproduce the results of the full 
form of the Bonn potential which has only one adjustable parameter: see Ref. [Hf 
for details. In the second, the OSBEP fll2f , mesons are treated as objects of non- 
linear theory. The mesons are the same as those in the Bonn B, but the form of 
the momentum space potentials is different. It is (we have taken into account that 
there is no energy transfer in the Hartree diagrams, i.e. u = 0) 

2 ~ (2pn + 1) 2 p"- 2 (1 + 2(1 - p)nf{S + \) n a^^ n 

v i{(l) — 9i ~ ~ ; ~ ~ 7> (5+l)(2pn+l) 



"=° {l + (4(p + 1R)-§(S + } 



1 



(24) 



q 2 + {2pn + l) 2 /i 2 ' 

where p = 1/2 for scalar mesons and p — 1 for pseudoscalar and vector ones, S is 
the spin of the meson, and the sum over n is practically converging at n = 4 [|l^ . 

Both these approaches permit one to check the status of the Walecka model |TJ 
by calculating the values of the vector and scalar fields in nuclear matter. For the 
case of charge-symmetric matter 

V = v w (0)p , S = - v al (0) + - A ^o(O)) Ps , (25) 
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where the scalar density p s is 



Ps 



p 



2t 



(2m + S- Vf 



g k FP 



(26) 



is the Fermi momentum and m is the free nucleon mass. Using the conventional 
equilibrium value of the nuclear matter density, p = 0.17 fm -3 , and the parameters 
of Table 5 of Ref.Q and Table 1 of Ref.@ we get 

V = +284 MeV , S = -367 MeV (27) 
for the Bonn B potential and 

V = +322 MeV , S = -404 MeV (28) 

for the OSBEP, both being close to those provided by the Dirac phenomenology flTSfl . 
So the contemporary iVW interaction potentials lead to nuclear relativity, the latter 
thus being really existing phenomenon rather than the suggestion of J.D. Walecka. 

For this reason the doorway state wave functions ij)\{x) should be treated as 
Dirac bispinors obeying the Dirac equation with 



n 



sp 



-«7o7V + i$(r)7 - + (to - l)m + V(r) + -? S(r) . (29) 
r 



The scalar and vector fields of finite nuclei consist of the isoscalar and isovector 
parts, the vector field also including the Coulomb potential 



S(r) = S (r) - T 3 Sx(r), r 3 
'3 



— 1, n 
+ 1, V 



, V(r) = V u {r) - r 3 V p (r) + ^^C(r), 



So(r) 
Si(r) 



v*i{q) + t v a0 (q) F s (q) e 



q 



(27T)* 



4 — - 4 
d 3 q 



d 3 r 



(27TY 



(30) 



C(r) 



Pch{r') 



r — r 



(27T) 



d 3 r' 



V P (r) 



v P (q)F (q)e 



(2tt) j 



$(r) - t — — 
2m dr 



where 



F s (q) = J Ps(r)e'^dh, F?(q) = J 



F(q) = J p(r)e-"^ 3 r, F~(q) = J ( 



P [r) + 



p;(r)e- iqr d 3 r 
k d 



{r 2 w-{r)) ] e~ 4qr d 3 r 



2mr 2 dr 

p{r) = p n {r) + p p {r), p~(r) = p n {r) - p p (r), p s {r) = p sn (r) + p sp (r) 
Ps( r ) = Psn{r) - p sp {r) , w~(r) = w n (r) - w p (r) . 



(31) 
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The scalar densities and the quantities w(r) are 



_ o(r) 2(T(r) + Hr) P >(r) + &(r)p(r)) 

Ps[r) ~ P[r) ~ (2m + S(r)-V(r))* [ } 

_ p>(r) + 2<S>(r)p(r) 

W{r) ~ 2m + S(r)-V(r) ' ^ 

They are calculated separately for neutrons and protons. The quantity r(r) is 
calculated in the local density approximation using Eq.(26). The isovector quantity 
<3>(r), Eqs.(29) and (30), arises from the tensor pNN coupling, k = f p /g p is the 
tensor-to- vector coupling ratio. 

So the two-particle contributions may be determined from experiment by using 
a definite model for the two-particle NN interaction. Little is known, however, 
about the many-particle NN forces. In such conditions it is reasonable to look for 
the many-particle contribution as a power series expansion over the nucleon density 
distribution: 

U m (r) = S m (r) + V m (r) = a 3 p 2 (r) + a 4 p 3 (r) + ••• (34) 

the p 2 (p 3 ) term resulting from three- (four)-particle forces etc. To elucidate the 
physical meaning of the coefficients let us consider a general form of the three- 
particle term: 

U 3 (r) = j h{\v-r 1 l\v-r 2 \)p{r l ) P {r 2 )d z v 1 d 3 r 2 . (35) 
In the homogeneous nuclear matter this gives 

U 3 = p 2 J h{\vv l l\v-v 2 \)dh l dh 2 , (36) 

and therefore 



In the same way 



as = J hMftfr, ■ (37) 

<H = //4(&^CWV S C- (38) 

These volume integrals are the only parameters which do not require any specific 
model for the many-particle NN forces. Such model is, however, necessary to take 
into account the finite range of the forces. We did not try to do this since (a) the 
problem of many-particle NN interaction mechanism is beyond the scope of our work 
and (b) the additional adjustable parameters describing the finite range cannot be 
safely determined because of the insufficient accuracy of the available experimental 
data. 
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The same reason forced us to introduce as little free parameters as possible 
and use all permissible simplifications. In particular, the many-particle terms are 
assumed to be equally distributed between the scalar and vector fields: 

S m (r) = V m (r) = \ U m (r) ■ (39) 

4 Results 

The observed and calculated spectra of the doorway state energies in Ca, 90 Zr and 
208 Pb nuclei are plotted in Figs. 4 and 5. 

The calculations are performed with two different two-particle potentials: the 
Bonn B, Fig.4, and the OSBEP, Fig.5. 

The results for the two-particle forces only are labelled as "pair". As seen from 
the figures the "pair" spectra are compressed compared to the observed ones, the 
lowest lsi/2 states being significantly underbound. This means that the potential 
well resulting from the two-particle forces only is too wide but insufficient deep, and 
so the actual well must be deeper and narrower as illustrated by Fig.6. 

Hence, the many-particle contribution (as discussed above this is the only reason 
for the difference between the actual and "pair" wells) consists of the attractive and 
repulsive parts, the radius of the former being less than that of the latter. The most 
simple form obeying this condition is provided by the sum of first two terms of the 
expansion (34) with a 3 > and < 0. In other words, the free-space many-particle 
NN interaction includes at least the three-particle repulsion and the four-particle 
attraction (of course the presence of higher many-particle forces is not excluded). 

Accounting for the fact that the many-particle forces contribute to both the 
isoscalar and isovector parts of the static nuclear field the quantity U m (r) is chosen 
in the form 

U m (r) = a 3 p 2 (r) + a 4 p 3 (r) - r 3 a 3 p(r) + a 4 p 2 (r) p'(r) 
p{r) = p n {r)+p p (r), p~(r) = p n {r) - p p (r) . (40) 

The finite size of nucleon is taken into account in the free-space NN forces, and 
therefore the static field of nucleus is expressed through the point nucleon densities. 
The proton ones p p {r) are obtained from the charge density distributions of Ref.j2|] 
by a usual deconvolution procedure. They are shown in Fig. 7a. The point neutron 
densities p„(r) are obtained from the folded densities of Ref. in the same way. 

The data of Ref. are based on high precision measurements of elastic electron- 
nucleus scattering thus providing the proton density distributions in the whole nu- 
clear region. The situation for the neutron densities is different since the elastic 1 
GeV proton-nucleus scattering underlying the data of Ref. || is sensitive mainly to 
the surface region of nucleus because of the absorption. For this reason the neutron 
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Figure 4: Spectra of doorway state energies with the Bonn B potential for the 
two-particle forces. 
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Figure 5: The same with the OSBEP. 
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U(r), MeV 

A 



Zr-90 




12 r, fm 



Figure 6: Isoscalar part of the static field in 90 Zr. The dashed and full curves are for 
the "pair" and actual wells respectively. The calculations are performed with the 
Bonn B two-particle forces and original nucleon density distributions of Ref. [3]. 
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Table 1: Neutron density parameters 





Bonn B 


OSBEP 


a 




a 


P 


40 Ca 


-0.0295 


0.5314 


-0.0255 


0.5230 


90 Zr 


-0.0758 


0.5551 


-0.0646 


0.5442 


208 pu 


-0.2645 


0.5445 


-0.2667 


0.5389 



density distributions p n {r) may differ from the Woods-Saxon-like ones of Ref. || in 
nuclear interior (as seen from Fig. 7a the proton densities are indeed different from 
the Woods-Saxon-like ones). The latter is just the region to which the doorway 
state energies are sensitive, and therefore they may be used to specify the Ref. || 
data for the neutron densities. We looked for the latter ones in the form 



Pn{r) = po 



W A (r) + aW A (0)<p 4 (M 



(41) 



where Wa{t) are the deconvo luted neutron densities of Ref. || and yj 4 (x) is the 
fourth Hermite function. The neutron density parameters a, (3 and the strength ones 
a 3; a 4 ; %; a 4 are determined from the best fit for both the doorway state energies 
and the elastic 1 GeV proton-nucleus scattering, the latter being calculated within 
the Glauber theory [|16|] . 

The density parameters are shown in Table 1. They are different for the two 
choices of the two-particle forces, but the difference is rather small. For this reason 
neither the resulting neutron density distributions, Fig. 7b, nor the 1 GeV proton- 
nucleus elastic scattering cross sections, Fig.8, are distinguishable in the figures. We 
also calculated the proton-nucleus cross sections with the original results of Ref. Q 
for the density distributions. As seen from Fig.8 the agreement with experiment is 
equally good for both the specified densities, Eq.(41), and the original ones. The 
many-particle strength parameters are 



5 

a7 = 25.5873 fm 5 , a7 = -128.5134 fm 8 (42) 



a 3 = 16.9296 fm° , a 4 = -107.6744 fm 

.5 



for the Bonn B two-particle forces and 

a 3 = 17.0011 fm 5 , a 4 = -110.3747 fm 8 , 

a7 = 26.9036 fm 5 , a7 = -130.1210 fm 8 (43) 

for the OSBEP ones. As seen from Eqs. (42) and (43) the strength parameters of 
the free-space many-particle forces are almost the same for the two cases. This is not 
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4\ 




2 4 6 8 10 12 U 16 18 0°C.m. 

Figure 8: 1 GeV proton-nucleus elastic scattering cross sections. The dashed and 
full curves are calculated with the original neutron density distributions of Ref. [3] 
and the specified ones respectively. 
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surprising because both the Bonn B and the OSBEP potentials provide an equally 
good description of the two-nucleon data, see the discussion in the Introduction. 

The results including both the contribution from the many-particle forces and 
the specified neutron densities are labelled as "full" in Figs. 4 and 5. The "full" 
doorway state energies agree with the observed ones (which are labelled as "exp") 
within the experimental error of 3 MeV. The exception is provided by the 2s\/2 
states in 208 Pb: in this case the discrepancy is about 5 MeV. The reason is not clear 
yet, but the discrepancy does not exceed two experimental errors. 

To estimate the relative importance of the two-particle and many-particle con- 
tributions to the static field of nucleus let us perform the calculations for nuclear 
matter, see Sect. 3. First consider the isoscalar part. The two-particle contribution 
is 



U 2 = V 2 + S 2 = he 



3 1 

4^1 (0) + |?Vo(0) ) ps 



-83 MeV, Bonn B 
-82 MeV, OSBEP , 
(44) 



whereas those from three-particle and four-particle forces are 



U 3 
U, 



hca^p 2 
hca^p 3 



96.5 MeV, Bonn B 
97 MeV, OSBEP 
-104 MeV, Bonn B 
-107 MeV, OSBEP 



(45) 



the many-particle contributions thus being as large as the two-particle one. 

The isovector part may be estimated by putting p~ = ^j^-p and p~ = 
The two-particle contribution is (see Sect. 3) 



Ps 



u; 



he 



v P (0)p~ - (v s {0) + \v al (0) - ^(0)) p s 



■N-Z 



MeV, Bonn B 
0.15^ MeV, OSBEP 



the many-particle one being 



(46) 



u: 



he (a 3 p + 



a 4 p 



(146 
(153-126 



125 = 21) ^ , Z MeV, 



27) ^ - Z MeV, 



Bonn B 
OSBEP 



(47) 



So the many-particle forces provide the dominant part of the isovector nuclear po- 
tential. The reason is due to the fact that the two-particle contribution arises from 
the exchange by isovector mesons p and S which are weakly coupled to nucleon, see 
Table 5 of Ref.|TH and Table 1 of Ref.lO. 
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5 Summary 



The above results give rise to the following general conclusions: 

1. Our results for the many-particle forces are quite competitive with those from 
the few-nucleon systems |L7| . Indeed, the properties of the latter ones (binding en- 
ergies, sizes, formfactors etc.) are expressed through the interaction in all orders of 
the perturbation theory, and therefore the solution of a rather complicated quantum 
mechanical problem is necessary to get the information on the many-particle forces. 
In contrast to the few-nucleon systems the doorway states for the one-nucleon trans- 
fer reactions in complex nuclei are solutions of a much more simple problem for one 
nucleon in central field. In addition the static nuclear field is expressed through 
the NN forces in first order of the perturbation theory, the results thus being very 
visual, see Figs. 1 and 6. The information from the doorway states is, however, re- 
stricted because it concerns only spin-independent terms of the many-particle forces 
(the spin-dependent ones do not contribute to the Hartree diagrams). Nevertheless 
it is a useful addition to that from the few-nucleon systems. 

Two important points should be mentioned in this connection, (i) Only three- 
particle forces (in addition to the two-particle ones) are included in all available 
calculations for the few-nucleon systems. Our results clearly show that this is in- 
sufficient, (ii) Calculating the nuclear correlation effects (binding energies and rms 
radii of finite nuclei, equation of state of nuclear matter, etc.), with the free-space 
NN interaction there is no reason to neglect the many-particle forces because they 
are as strong as the two-particle ones, compare Eqs. (44) and (45). 

2. The effective three-particle and four-particle forces are also repulsive and 
attractive respectively in the recent calculations within the relativistic mean-field 
approximation [|18|, 0, see the Appendix. Such forces include implicity the cor- 
relation effects which are not taken into account explicitly within this framework. 
For this reason the above signs of the forces might be treated as the artifact of the 
approximation. But our results for the free-space many-particle forces show that 
this is not the artifact. 



A Appendix 



The potential energy of the a mesons is fllS , |I~9 



v ' 2 3 4 

with A 3 < and A 4 < 0, the scalar field S = ga thus obeying the equation 



A 



3 „3 



A 



4 -4 



A, 



A 



(A.l) 



(A.2) 
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Let us use the following iteration procedure 

(A - » 2 )S n = gl Ps + + (A.3) 

9 9 

with 

(A-^)S = glp s (A.4) 
for the initial iteration. The result is 

S(r) = -g 2 jydv-r^i^dh, + 

+ / /3(|r - ri|, |r - r 2 \)p s (n)p s (r 2 )dr 1 dr2 + (A. 5) 

where 

= (A.6) 
47r|r — r'| 

/ 3 (|r- ri |,|r-r 2 |) = -A 3 £ 3 / y(|r - r'|)y(| ri - r'|)y(|r 2 - r'|)dV (A.7) 
/ 4 (|r - n|, |r - r 2 |, |r - r 3 |) = 
= A^/ydr-r'IMIn-r'IMI^-r'Dydra-r'DdV - (A.8) 

- 2X^1 y(\r - r'IMIn - r'\)y(\v' - r"\)y(\v 2 - r"\)y(\r 3 - r"|)Wr". 

The dots in the rhs of Eq.(A.5) represent the higher-power terms in respect of p s 
resulting from the higher many-particle forces. As seen from (A.7) the three-particle 
force is repulsive because of the sign of A3 (g > in Refs.[|l8], |1|). The four-particle 
one, Eq.(A.8), consists of two terms. The first is of first order with respect to the 
A4 term of (A.l). It is attractive because of the sign of A4. The second is of second 
order with respect to the A3 term. It is attractive irrespective of the sign of A3. 
The volume integrals of the forces (A.7) and (A.8), Eqs. (37) and (38), are 

«s = -2 3 A 3 /i- 6 , a 4 = # 4 (A 4 - 2X 2 3 ^ 2 )p~ s . (A.9) 

The least values of these quantities correspond to the NL-SH parameter set of 



Ref.|T8[, Table 2 of this reference. They are 

a 3 = 21.9 fm 5 , a 4 = -136.5 fm 8 , (A.10) 
thus being rather close to the free-space values, Eqs. (42) and (43). 



20 



References 

[1] M. Baranger, Nucl.Phys. A149 (1970) 225. 

[2] H. de Vries et al. At. Data and Nucl. Data Tables 36 (1987) 495. 
[3] G.D. Alkhazov et al. Nucl.Phys. A381 (1982) 430. 

[4] A. A. Abrikosov, L.P. Gor'kov and I.E. Dzyaloshinski "Methods of quantum 
field theory in statistical physics" M, 1962. 

[5] A.B. Migdal "Finite Fermi-system theory and properties of atomic nuclei" M. 
1983. 

[6] S.G. Kadmenski and PA. Lukyanovich, J.Nucl.Phys. 49 (1989) 1295. 
[7] S.S. Volkov et al., J.Nucl.Phys. 53 (1990) 1339. 
[8] A.A. Vorobyov et al., J.Nucl.Phys. 58 (1995) 1923. 
[9] K. Erkelenz, Phys.Reports 13 (1974) 191. 
[10] M. Lacombe, Phys.Rev. C21 (1980) 861. 

[11] R. Machleidt, K. Holinde and Ch. Elster, Phys.Reports 149 (1987) 1. 

[12] L. Jade and H.V. von Geramb, Phys.Rev. C57 (1998) 496. 

[13] R. Machleidt, Adv. in Nucl.Phys. 19 (1989) 189. 

[14] J.D. Walecka, Ann. of Phys. (NY) 83 (1974) 491. 

[15] S.J. Wallace, Comments on Nucl. and Part. Phys. 13 (1984) 27. 

[16] R.J. Glauber, "Lectures in theoretical physics", ed. W.E. Brittin et al., vol.1 
(NY 1959). 

[17] R. Schiavilla, V.R. Pandharipande and R.B. Wiringa, Nucl.Phys. A449 (1986) 
219. 

[18] G.A. Lalazissis, J. Kdnig and P. Ring, Phys.Rev. C55 (1997) 540. 
[19] M.L. Cescato and P. Ring, Phys.Rev. C57 (1998) 134. 



21 



